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The growth of N,N'-dimethylperylene-3,4,9,10-bis(dicarboximide) (DiMe-PTCDI) on KBr(OOl) 
and NaCl(OOl) surfaces has been studied. Experimental results have been achieved using frequency 
modulation atomic force microscopy at room temperature under ultra-high vacuum conditions. On 
both substrates, DiMe-PTCDI forms molecular wires with a width of 10 nm, typically, and a length 
of up to 600 nm at low coverages. All wires grow along the [110] direction (or [110] direction, 
respectively) of the alkali halide (001) substrates. There is no wetting layer of molecules: Atomic 
resolution of the substrates can be achieved between the wires. The wires are mobile on KBr surface 
but substantially more stable on NaCl. A p(2 x 2) superstructure in brickwall arrangement on the 
ionic crystal surfaces is proposed based on electrostatic considerations. Calculations and Monte- 
Carlo simulations using empirical potentials reveal possible growth mechanisms for molecules within 
the first layer for both substrates, also showing a significantly higher binding energy for NaCl(OOl). 
For KBr, the p(2 x 2) superstructure is confirmed by the simulations, for NaCl, a less dense, 
incommensurate superstructure is predicted. 

PACS numbers: 68.37.Ps, 61.64,+w, 68.55.Ac 



I. INTRODUCTION 

The possible use of organic molecules for advanced 
types of electronic devices has attracted much interest 
in recent decades. [TJ [2j [3] A large field of research has 
emerged, studying the structural and electronic prop- 
erties of various molecules on solid surfaces. Scanning 
probe methods have proved to be a powerful tool for these 
studies, especially scanning tunneling microscopy (STM) 
has been used for many promising studies. [4] [5j [6] Nev- 
ertheless, STM is limited to conducting surfaces, while 
insulating surfaces are still important: Here, the influ- 
ence of the substrate on the electronic properties is small 
compared to metal or semiconductor surfaces. 

Frequency modulation atomic force microscopy (FM- 
AFM)[7J |8] is a scanning probe technique that allows 
for real space imaging of insulating surfaces. Thus it 
has made the analysis of the structural properties of 
organic molecules on insulators accessible. [9j [10] Re- 
cent studies demonstrate the capability of the tech- 
nique to analyse the growth of the prototypical semi- 
conductor molecules C6o[ffl D2] and 3,4,9, 10-perylene- 
tetracarboxylic-dianhydride (PTCDA) [El HI [111 M HZ] 
on alkali halide (001) surfaces. 

However, all these studies show dewetting of the 
surfaces and the growth of bulk-like molecular struc- 
tures for PTCDA or fractal islands for Ceo- For 
the development of small-scale electronic devices, a 
controlled growth of low-dimensional molecular struc- 
tures would be desirable. [18] This kind of growth 
has recently been analysed for a perylene derivate 
similar to PTCDA, N,N'-dimethylperylene-3,4,9,10- 
bis(dicarboximide) (DiMe-PTCDI); [EJ EQJ EH [22] this 
molecule shows the growth of wire-like structures on 



* email: markus.fendrich@uni-due.de 



KBr(001).[23] 

The present paper reviews the results for KBr (001) and 
shows experimental data for another alkali halide surface, 
NaCl(OOl), as well as molecular force field calculations 
analysing the molecular arrangement of DiMe-PTCDI on 
both substrates. 



II. EXPERIMENTAL 

The experiments have been performed at room tem- 
perature under ultra-high vacuum(UHV) conditions in a 
commercially available analysis system (AFM/STM by 
Omicron, Germany) operated in FM-AFM mode. The 
oscillation of the cantilever is detected by beam deflec- 
tion; to illuminate the cantilever, a SuperLum LED 
coupled to the UHV chamber by a single-mode fiber 
was used. For the detection of the frequency and con- 
trol of the cantilever oscillation amplitude, we used the 
easyPLLplus system (Nanosurf, Switzerland), operating 
in self oscillation mode. The specifications of the sili- 
con cantilevers (PPP-NCL by Nanosensors, Switzerland) 
were: resonant frequencies of fo ~ 160 kHz, spring con- 
stants of k « 30 - 40 N/m and quality factors of Q « 
30000 in UHV. The cantilevers were driven at an am- 
plitude of A = 20 nm, typically. In some experiments, 
a voltage between tip and sample was applied to com- 
pensate for electrostatic interactions. Specific frequency 
shifts are given for each image. 

For scanning control and data acquistion, the open- 
source software GXSM[24 combined with home-built 
electronics has been used. Image processing has been 
done using the free software WSXM.[25] 

Clean KBr and NaCl (001) surfaces have been prepared 
by in situ cleaving of commercially available single crys- 
tals (Korth, Germany). The crystals have been heated at 
450 K for 1 h before cleaving to remove contaminations 
and also after cleaving to remove trapped charges. After 
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preparation of the substrates, various amounts of DiMe- 
PTCDI have been evaporated onto the samples from a 
home-built crucible. The temperature of the oven was 
550 K, the evaporation rate was checked before and after 
deposition by a quadrupole mass spectrometer. During 
evaporation, the sample was kept at room temperature. 



III. EXPERIMENTAL RESULTS 
A. KBr(OOl) 

On KBr(OOl) surfaces, DiMe-PTCDI grows wire-like 
islands with a length of up to 600 nm (see Fig. [TJ. Be- 
tween the wires, atomic resolution of the KBr substrate 
could be achieved: There is no wetting layer of organic 
molecules. All wires are aligned with (110) surface axes. 
The wires are stabilized by step edges, "loose ends" are 
seldom observed. When no step edges are present, the 
structures are highly unstable and diffuse when scanned 
for a longer time. The wires have a height of at least two 
molecular layers. 

Due to their instability, no molecular resolution of the 
wires was possible, so no information on the arrangement 
of the molecules within the wire, could be achieved exper- 
imentally. For the initial growth, i.e., for the first layer 
upon the KBr substrate, a model proposing a p(2x2) 
superstructure in brickwall arrangement has been elabo- 
rated. This model was based on force field calculations 
and electrostatic considerations for a single molecule on 
the ionic KBr(OOl) surface. [23] 




FIG. 1: (color online) [25] Review of the previously reported 
experimental results for DiMe-PTCDI on KBr(001).[231 The 
organic molecules grows wire-like structures with height of at 
least 2 molecular layers at low (a) and higher (b) coverages, 
all wires aligned with (110) surface axes. When the sample is 
cooled during evaporation (c), the height of the wires is only 
one layer. The wires are stabilized by substrate step edges; 
when no step edges are present, the wires are unstable and 
diffuse, as shown in the series of images (d)-(g). All images 
have been obtained at small frequency shifts, df = -0.5 Hz., 
image sizes are (a) 600 x 600 nm 2 , (b) and (c) 400 x 400 
nm 2 , (d)-(g) 300 x 300 nm 2 . 



B. NaCl(OOl) 

On NaCl(OOl) surfaces, DiMe-PTCDI molecules show 
similar structures. As shown in Fig. [2j wires are still 
aligned with (110) surface axes. No wetting layer of 
molecules is found. In contrast to KBr, the wires are 
much more stable, even if they are not stabilized by step 
edges: The area shown in the image has been scanned 
for several hours without showing any changes in the ar- 
rangement of the molecular islands. All the wires are two 
molecular layers high. 

On KBr, continuous networks of wires have been 
found, supposably due to stability reasons; here, we find 
smaller islands that are not interconnected and some- 
times consist of several merged perpendicular wires. No 
step edges are needed as nucleation centers for the is- 
lands on NaCl; it seems reasonable that islands nucleate 
at surface defects. 




FIG. 2: (color online) [25] Molecular wires of DiME-PTCDI 
on NaCl(OOl) at low coverage (0.2 ML). The image has been 
obtained at a frequency shift of df = -0.5 Hz, the size of the 
image is 400 x 400 nm 2 . No NaCl step edges can be seen 
within the image. All wires are oriented towards (110) surface 
directions and very stable. No changes could be observed 
while scanning this area for several hours. There is no wetting 
layer of molecules, between the wires, atomic resolution of the 
substrate is possible (see inset, 3x3 nm 2 , df = -6 Hz). 

When the coverage of molecules on the surface is in- 
creased, larger and more compact molecular islands are 
found while the framework of the perpendicularly ar- 
ranged wires can still be seen. Islands with a height of up 
to 4 molecular layers have been observed. Within these 
islands, 1-layer areas are embedded; however, these wet- 
ting layers are unstable, a dewetting of the 1-layer ar- 
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eas is observed within several minutes. Fig. [3] shows 
two subsequent scans of the same area (500 x 500 nm 2 ). 
One can see the molecules dissappear from the enclosed 
areas covered with only one layer of molecules. These 
molecules diffuse two higher layers, as reported also for 
C 60 molecules on KBr(001).[Il] H2] 

As the islands are wider and more stable than on 
KBr(OOl), molecular resolution of the topmost layer of 
DiMe-PTCDI on NaCl(OOl) was possible. Figg shows 
resolution of molecular rows on top of an island which is 
2 layers high and merged out of 3 wires. The junction 
can clearly be seen on the molecular scale as an inter- 
face of perpendicular rows. The internal structure of the 
molecular rows cannot be resolved; nevertheless, every 
second row seems shorter at the junction, thus a brick- 
wall arrangement of the molecules seems reasonable. 




FIG. 3: (color online) [25] Images (a) and (b) show two sub- 
sequent scans of a 500 x 500 nm 2 large area at high coverage 
(the total coverage is about 2 ML of DiMe-PTCDI). The time 
between (a) and (b) is about 30 min. One can see large is- 
lands of molecules with heights between 1 and 4 molecular 
layers. The rectangular structure aligned with (110) surface 
axes can still be seen, in contrast to the situation at low cover- 
ages, the islands are more compact and do not show wire-like 
growth. The ML-covered areas (see arrows) are not stable, 
the dewetting of these areas and some changes in the shapes 
of the higher areas can be seen in (b) . Both images have been 
obtained at a frequency shift of df = -0.5 Hz. 



IV. CALCULATIONS 

To get further insight into the growth mechanisms for 
DiMe-PTCDI on alkali halide (001) surfaces and corro- 
barate the models based on the experimental findings, 
molecular force field calculations based on empirical po- 
tentials have been performed. The capability of em- 
pirical potentials has been shown for several molecular 
systems. [26, 27, 28, 29 In particular, similar calculations 
could elucidate the experimental data for a very similar 
system, PTCDA on KCl(001).[l6] 

As preliminary calculations, the potential energies 
of single DiMe-PTCDI molecule on KBr(OOl) and 
NaCl(OOl) have been calculated using the empirical po- 
tentials of the AMBER molecular force field. [13 EEJ [32] 



FIG. 4: (color online) 25 The left image (100 x 100 nm 2 , 
df = -0.5 Hz) shows a 2-layer high island of DiMe-PTCDI on 
NaCl(OOl). The island is merged out of several perpendicular 
wires, the right image shows a zoom (12 x 12 nm 2 , df = - 
5.5 Hz) of the junction. Molecular rows can be resolved; at 
the junction, every second row appears shorter (see arrows). 
Therefore, a brickwall arrangement of the molecules seems 
reasonable. 



As DiMe-PTCDI is a polar molecule (the electronega- 
tivity of the 4 oxygen atoms causes a shift of negative 
charge to the "corners" of the molecule), the electrostat- 
ics of the system have to be accounted for. The polarity 
of the molecule is represented by partial charges assigned 
to each atom within the molecule by a Mulliken method, 
using the software Hyperchem. Each ion of the substrate 
is represented by a point charge of ±eo- 

Calculating the potential energy as described above, 
the minimum energy position for a single molecule was 
determined for both substrates using a gradient search 
algorithm. Optimum positions are found as shown in 
Fig. [5j i.e., on top of a halide ion, the molecule aligned 
with the [110] axis of the substrate. In this way, the 
negatively charged oxygen atoms are close to 4 positive 
alkali ions of the substrate. The adsorption geometry 
is more favorable for NaCl, which has a smaller lattice 
constant of 562 pm compared to KBr (658 pm); therefore, 
the binding energy of a single molecule is higher for NaCl 
(E NaCl = - 1.65 eV, E KBr = -1.35 eV). 

In a second step, potential energies for larger molec- 
ular domains have been calculated. The Grid method 
and the software Powergrid [26j [33] have been used 
to calculate the potential energy for domains consist- 
ing of 150 molecules or more. The molecular arrange- 
ment within the domain has been optimized using a 
Monte-Carlo method, the Metropolis algorithm. [34 Sev- 
eral Monte-Carlo experiments have been performed: Af- 
ter the Metropolis algorithm was started for an arbi- 
trary initial configuration, the arrangement was opti- 
mized, i.e. the size of the unit cell and the orientation of 
the molecules within the unit cell were changed until an 
optimum arrangement was found. The calculation con- 
firms the supposed p(2 x 2) brickwall arrangement for 
one layer of DiMe-PTCDI molecules on KBr(OOl). For 
NaCl(OOl), a incommensurate structure with an epitaxy 
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FIG. 5: (color online) Minimum energy positions for DiMe- 
PTCDI molecules on KBr(OOl) and NaCl(OOl), calculated 
with an AMBER force field calculation. In both cases, the 
minimum energy position is aligned with the [110] substrate 
axis, on top of a halide ion (Br - or CI - ). Due to the smaller 
lattice constant of NaCl, the binding energy is higher for this 
substrate(i?jVaCJ = -1.65 eV, Ekbv — -1.35 eV). 



matrix of 
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was found. Within the domain, all the molecules are still 
aligned towards the (110) substrate axis. 

The higher stability of the islands on NaCl(OOl) reflects 
in a higher binding energy for the domains. The energy 
per molecule within a domain is Ekbv dom — ~ 

1.4 eV for 

KBr and E Na ci,dom = -1-8 eV for NaCl. Table [I] shows 
the results of the calculations. 

TABLE I: Summary of the results of the Monte Carlo simu- 
lations for DiMe-PTCDI on KBr(OOl) and NaCl(OOl). 





KBr(OOl) 


NaCl(OOl) 


structure 


commensurate 


incommensurate 




p(2 x 2) 






brickwall 


brickwall 


area per 


A = 1.74 nm 2 


A = 1.37 nm 2 


molecule 






binding energy 


-1.35 eV (single) 


-1.65 eV (single) 


per molecule 


-1.4 eV (domain) 


-1.8 eV (domain) 



V. CONCLUSION 

We have shown an experimental study of DiMe-PTCDI 
on two alkali halide (001) substrates, KBr (001) and 
NaCl(OOl). On both substrates, we have observed the 
growth of wire-like molecular islands which are oriented 
along the (110) substrate axis. The wires are at least 
2 molecular layers high for both substrates. The wires 
are very mobile on KBr surfaces, where they are sta- 
bilized by step edges; on NaCl surfaces, the wires are 
more stable, allowing for the resolution of molecular 
rows within the topmost molecular layer. When larger 
amounts of molecules are evaporated onto NaCl(OOl) sur- 
faces, large, compact islands are found. In contrast to the 
two- layer wires at low coverages, monolayer- covered ar- 
eas are found. However, a dewetting of these areas was 
observed, possibly induced or supported by the scanning 
tip. 

Calculations based on empirical potentials give further 
insight to the structural properties. For KBr (001), a p(2 
x 2) superstructure is found, all molecules are aligned 
towards the (110) substrate axis. This alignment is also 
found for NaCl(OOl), however, in this case an incommen- 
surate superstructure is predicted. However, at least this 
purely theoretical finding has to be doubted, as there are 
no experimental hints. Inaccuracies of the calculation 
may occur for several reason: In contrast to the exper- 
iment, only one layer of molecules on the substrates is 
considered; all molecules are lying flat on the surface, a 
constriction that is necessary for the Grid method; [26] 
and possible conformational changes of molecules or sub- 
strates are completely neglected, both molecules and sub- 
strates are rigid. In spite of the inaccuracies of the calcu- 
lations, the higher stability for NaCl reflects in a higher 
binding energy for this substrate compared to KBr. 
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